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Abstract 

The deuteron form factors are calculated in the framework of the relativistic nucleon- 
^3 ! meson dynamics, by means of the explicitly covariant light-front approach. The inflluence 

O^ ' of the nucleon electromagnetic form factors is discussed. At Q^ < 3 (GeV/c)^ the predic- 



a^ 



tion for the structure function A{Q'^) and for the tensor polarization observable t2o are in 
agreement with the recent data of CEBAF/TJNAF. 



§ : 1 Introduction 

- Recently the first experimental data obtained at Thomas Jefferson National Accelerator Facility 

r^ I (former CEBAF) on the deuteron structure function A and tensor polarization ^20, measured 

c^ I in the elastic ed scattering, have been reported |]T|-Q. They are rather precise and correspond 

to a maximum momentum transfer of 6 (GeV/c)^. Since in the ed collision the deuteron gets 

from the electron a momentum comparable with the deuteron mass, these data are probing the 

truUy relativistic dynamics inside the deuteron. 

Our preliminary calculations of the form factors and of ^20 have been presented in [|]. 
We present in this paper the details of our approach, discuss the influence of the different 
choices of nucleon electromagnetic form factors parametrization and compare our results with 
the reported experimental data. 

In our calculations we assume the deuteron structure to be determined by the relativistic 
nucleon-meson dynamics. Namely, we suppose that the nucleons in the deuteron interact by 
exchanging relativistic mesons. We take the same set of mesons and parameter values used 
in the construction of the Bonn potential |0. However, we do not make any nonrelativistic 
potential approximation and calculate the one-boson-exchange kernel in relativistic form, as it 
appears from field theory. 

In the momentum tranfer region scanned in TJNAF, the relativistic effects related to 
the nucleons motion, to the spin rotations and to the retardation of the exchanged mesons, 
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should manifest themselves in full measure and be of crucial importance in describing the 
data. Therefore we believe that their main part can be properly taken into account even in an 
approximate relativistic calculation. On the other hand, it is important to work in an approach 
which provides a clear physical interpretation of the incorporated effects. 

The deuteron form factors are calculated in the framework of the explicitly covariant 
version of the Light-Front Dynamics (LFD) recently reviewed in ^j. In this approach, the state 
vector is defined on the light-front plane of general position uj-x = 0, where u; is a four-vector 
such as u"^ = 0. This restores the relativistic covariance lost in the standard light-front approach 
which is obtained as a particular case for uj = (1,0,0, —1). The relativistic wave functions - 
the Fock components of the state vector - are the closest couterparts of the nonrelativistic 
ones. This allows one to benefit from the knowledge of the nonrelativistic wave functions 
and to incorporate selfconsistently the relativistic effects. Using the light-front plane leads to 
significant simplifications due to the absence of vacuum fluctuations. The diagram technique 
is a three-dimensional one and can be interpreted in terms of time-ordered physical processes. 
Therefore the selection of the diagrams contributing to the kernel or to the form factors, whose 
numerical estimation is rather difficult in any approach, can be supported by the intuitive 
understanding of the corresponding physical process. All that allows, in a given relativistic 
dynamics, to carry out the calculations on a satisfactory level of confidence. 

We mention also few other approaches to calculate the relativistic wave function and the 
deuteron electromagnetic form factors. In the standard version of LFD, defined on the plane 
t + z = 0, these form factors were calculated in ||7|-p^. The results obtained in |11, 12] are based 



on the Bethe-Salpeter approach and the works |]13l0] use a three-dimensional reduction of the 



Bethe-Salpeter function - the Gross wave function [|T^ . 

In sect. |] we present the relativistic deuteron wave function on the light front. In 
sect. ^ the electromagnetic vertex of the deuteron in the impulse approximation, based on 
the one-body electromagnetic nucleon current, is discussed. In sect, ^we discuss the contact 
(instantaneous) interaction, corresponding to the NNB'-f vertex {B is an exchanged meson) 
which is a correction to the one-body current. In sect. |^ we explain how to extract the 
form factors, separating them from the nonphysical contributions. Section ^ is devoted to the 
evaluation of the nucleon electromagnetic form factors infiuence in the deuteron observables. 
Section ^ contains a comparison of our results with the experimental data and concluding 
remarks. 

2 Wave function 

The wave functions are the Fock components of the state vector defined on the light-front plane 
uj-x = 0. The explicit covariance allows one to construct the general form of the light-front 
wave function for a system with a given spin. The relativistic deuteron wave function on the 
light front contains six spin components, in contrast to two components - S and D-waves - 
in the nonrelativistic case. Its general form is given in |TB|,|r^. Below we will keep only three 
dominating components: 



!Z/f-i = ,/^e'i{p)r^ik2)(j)McU'''ik 



1^ 



{ki -k2)u ,1 i ;a;/3 7 /i ^ 

/m^ m m'^uj-p 



where p and ki^2 are the on mass shell deuteron and the nucleon momenta, e^(p) is the deuteron 
polarization vector, u'^{k) is the nucleon spinor, Uc is the charge conjugation matrix. We notice 
that the wave function defined on the light-front plane depends on the orientation of this plane 
through the argument u. 

In the system of reference where /ci + /c2 = the function (|I|) obtains the more transparent 
form: 

Mn) = h^^ + f^^^-^ - -) + f^\[llf X nl (2) 

— * — * 

where k is the value of ki in this system of reference, n is the direction of cu in this system, Wg. 
is the two-component nucleon spinor. The scalar functions fi depend on the scalars k = \k\ and 
z = n-k/k. One can find a special representation [jl^ in which the wave function (|l]) obtains the 



form (|^) in arbitrary system of reference. The components yjj in (|I|) and fi in (0) are linearly 
related with each other P,p!7|. 



The wave function (|T]) was calculated in [|T7| in a perturbative way incorporating the full 
relativistic one-boson exchange (OBE) kernel found in LFD, with the corresponding nonrela- 
tivistic wave function as zero order approximation. This calculation can be considered as the 
perturbation theory developed in terms of the difference between the relativistic kernel and 
the nonrelativistic potential. For the OBE kernel the set of mesons, coupling constants and 
form factors corresponding to the Bonn model P were used. The solution thus obtained is 
approximate in two ways: first for its perturbative character does not provide an exact solution 
of the LFD equation and second because the parameters of the OBE kernel were kept to their 
original values i.e. were not fitted in the LFD relativistic framework. The accuracy of our 
solution in reproducing the low energy deuteron observables was estimated at the level of 20% 
in norm. 

In nonrelativistic region of k, components /i and /2 found this way turn into the usual 
S- and D-waves, whereas other components become negligible. However, starting from k p^ 0.5 
GeV/c, component /s dominates over all other components, including /i and /2. 

The physical meaning of this dominating extra component has been clarified [|18| by 
comparing, in l/m approximation, the analytical expression for the amplitude of the deuteron 
electro disintegration near threshold with the nonrelativistic one, including meson exchange 
currents. For the isovector transition, in the region where the so called pair term with the 
pion exchange dominates, this component (together with a similar component in the scattering 
state) automatically incorporates 50% of the pair term contribution and therefore dominates 
too. Another 50% is given by the contact interaction (see below). 

3 Electromagnetic vertex 

The electromagnetic vertex J^^, we use to calculate the deuteron form factors corresponds to 
the sum of the impulse approximation (lA) (figure |l|) and of the contact interaction (C) (figure 
^ which takes partially into account the two-body current: 

J,^, = J,^,(IA) + j;,(C). (3) 

These contributions are explained in what follows. 



The deuteron lA electromagnetic vertex is shown in figure ^ The dashed fine is associated 
with a fictitious particle - the so called spurion. It refiects the fact that, although all the 
momenta are on the corresponding mass shells, the wave function is off energy shell, and 
therefore there is no any conservation law between all the components of the nucleon and the 
deuteron four- momenta. The spurion momentum just absorbs their nonzero difference and one 
has ki + k2 — p = OJT. To avoid misunderstanding, we emphasize that the spurion line does not 
imply the presence of an extra particle in the intermediate state in figure |I| which contains two 
nucleons only. 




Figure 1: Electromagnetic vertex of the deuteron. 



The electromagnetic amplitude in the impulse approximation is derived by applying the 
rules of the graph techniques ||^ to figure |T]. It has the form: 

(AV''(IA)|A)=eI^(p')J,^.(IA)e^(p), (4) 

Jjlum = J^sJ TrW, {k', + m) r {k, + m) Mh - m)] ^^ _ ^y^^^^ , 

where x = ui-ki/uj-p. The wave function 0,^ = (f)i^{ki,k2,p,ujT) given by eq.(|l]) corresponds to 
the deuteron initial state, while 0^ = 0^(fci, k2,p',ujT') corresponds to its final state. Tp is the 
electromagnetic vertex of the nucleon: 

r^ = Fi7^ + ^a''-g„, (5) 

2m 

Fi and F2 being the nucleon elecromagnetic form factors. 

The expression for the trace in (^ is obtained as follows. To each nucleon line in figure |l| 

we associate the LFD nucleon propagator, with the NN'j vertex we associate expression (j^), 

and with d — NN vertices the deuteron wave functions. We get in this way the following product 

of spin matrices: 

e^ip') {bo^'.UcloHK + m)TPCk2 + m)4>Mc}^^ e\{p){k^ + m)^„. (6) 

Factor {. . .j^^ in (H) corresponds to the upper line of the diagram and the factor (A;i+m)^Q, 
corresponds to the lower line. We attach the deuteron wave function to the upper one. We 
keep in (P) the matrix indices a and j3 explicitly. Since both nucleon lines are passed in the 



same direction, the order of indices /3, a is the same. This means that one of the factors in (|^) 
(we take the second one) is a transposed matrix. The factor [7o0uf^c7o]^ originates from the 
conjugated final deuteron wave function: 



where = 7o0^7o- With the wave function (|I|) we get: 

bo4>',uM^ = -t/,0;. 



(7) 



and, hence, obtain the factor: 



-Uc{ki + mfUc = {ki - m). 



that gives the trace in eq.(^). A more detailed derivation of (^) by means of the LFD graph 
technique can be found in ^j . Next section is devoted to the derivation of the contact contri- 
butions. 



4 The contact interaction 

A peculiarity of the LFD is the existence, in addition to the impulse approximation, of the so 
called instantaneous (or contact) interaction in the NNB'j vertex, where B is any of the mesons 
7r,p, . . ., building the NN potential. Its contribution to the electromagnetic vertex in g^ order 
is shown in figure |^. The cross on a fermion line means that this line is not to be associated 
with a propagator, but with a factor proportional to a; = a'^7'^. In the standard approach this 
is the well known instantaneous interaction vertex 7^ The nucleon electromagnetic vertex Tp 
is determined, as usual, by the two nucleon electromagnetic form factors given in eq.(||). 






(b) 




Figure 2: Contact term contributions to the electromagnetic interactions with deuteron. 



In evaluating the contact amplitude, we have taken into account the sum over all six 
mesons contributing to the Bonn potential, with the parameters used in the Bonn model [^]. 



This amplitude incorporates partially the two-body currents. In principle more complicated 
diagrams, of order higher than g"^, may also contribute. Since our wave function takes into 
account corrections associated with the first degree of (7^, we keep the same order for the 
electromagnetic vertex. This means that we consider only the diagrams of figure |] and in 
calculating these diagrams we omit component /s in the wave function, since its contribution 
multiplied by the contact interaction is of higher order than g^. To calculate the effects beyond 
the g^ order, one has also to take into account the higher order irreducible contributions to 
the OBE kernel. One can expect that these contributions are incorporated in this kernel on 
a phenomenological level by an appropriate fit of its parameters (coupling constants, meson 
masses and cutoff parameters). However, this is not the case for the electromagnetic vertex. 

The contact amplitude is given by the sum of two terms which differ by the relative time 
order of the contact and the electromagnetic vertices: 

JU^) = J',AC)\left + J',AC)U,M. (8) 

The Jl^ui^)\ieft term corresponds to the sum of diagrams (a) and (b) in figure |, which 
differ from each other by the relative time order of the meson vertices. Both of them have the 
contact vertex on the left hand side of the electromagnetic one. The opposite order is assumed 
in Jl^ui^)\right which corresponds to diagrams (c) and (d) in figure ^. In addition one has to 
take the sum over all mesons contributing to the interaction. 

For a given meson, the amplitude corresponding to the sum of diagrams (a) and (b) in 
figure - left contact term - has the form: 



Jpieft = mJTT{{-U,4>',){k'2 + m)TP 




^^ ' ^ ^ '^\^'' + K^) (27r)3(l - x)2ek, (27r)3(l - x')2ek'^ ' ^ > 

where x = oo-ki/uo-p, x' = uj-k[/uj-p. The factor (/i^ + K"^) is the denominator of the meson 
propagator, which is expressed in terms of the relative nucleon momenta. In these variables K^ 
is given by 



K^ = ik' _ kf - in-k'){n-k) ^''' "'^' + (4 + 4 _ Im^) 



^k'^k 



n-k ' n-k 



^k' ^k 



(10) 



In the time ordered graph technique, the analytical expressions for the diagrams (a) and 
(b) in figure ^ differ by the meson propagators. The modulus in (p!OD just take into account this 
change of the meson propagator for different time ordering of the vertices. Therefore expression 
(y) with (|T0|) for K'^ corresponds to the sum of graphs (a) and (b). 

Expression (^ is obtained similarly to the case of the impulse approximation. The factor 
V stands for the meson-nucleon vertex at the upper and lower lines. The factor {—Uc4>'n) 
appears from relation (|^). Taking again into account the fact that Uck^ = —kUc we get: 

Tfl r ^ A A -N 

'^'^■'^''^' = 2H^iuj-p) J ^""W^'^ + ^)^'^V{k, + m)ct>,{k, - m)V\k[ - m)} 

1 1 d^kcPk' 



{li^ + K^){l-x') Sk Ek 



where V^ = UcV^Uc- We used the relation: uj-k2 = (1 — x')uj-p. Instead of momenta ki, k2 we 
integrate in (|TT|) over the relative nucleon momenta k, k ', using the relation: 



d^ki 



(fk cPR I dx 



2(l-a;)efci £k 2s(l - a;) 



(12) 



and similarly for k[,k', R'^. These relative momenta are the arguments of the initial and final 
deuteron wave functions respectively. 

The right contact term is obtained from eq. (|ll|) by the replacement 



TpQjV 



VujTP 



(1-x') (1-x) 



(13) 



In case of scalar (pseudoscalar) exchanges one should put both for up and down vertices 
V = V^ = g {y = V^ = ig'y^)- For the pseudoscalar exchange one can simplify the contact 
terms by excluding 75: 



jp \ps 



m 



277r6(cu-p) 



Tr{(f)'^{k2 + m)T'^Lj{k2 — 'm)(f)^{ki + m){k'i — m)} 

g2 1 d^k d^k' 



{fl^+K^)il-x') Ek Sk' 



(14) 



JP |P« 

fiuiright 



m 



2^6 (^.p) 



Tr{(j)' (k'2 — m)ujV^{k2 + m)(j)y(ki — m)(k[ + m)} 



d'^k d^k' 



X 



(yU2 + i?2)(l-a;) Ek Sk' 



(15) 



where (f) = 75(/'75. The function (j) differs from (j), eq.(0), by changing the sign of (^2 (and also of 
(Pi and (Pq if they are not neglected) . 

Let us now find the vertices V and V^ for the vector exchange. It is convenient to take 
them from the expression of the corresponding kernel which has the form [^ : 



K 



^2(91" - -— o-"'"z(A; - ujri)a')u2 
2m 



-gaf3 + 



{k - UJTi)a{k - UJTi)f3 



fX^ 



2m 

dr 

x5{{ki -k[ + uTi - urf - n^)e{ujiki - k[)) — 

Ti — le 

{k-UJTi)a{k-UJTi)p 



+ ju'^igT + ^^"'"^(^ - ^n)a')w2 



-9al3 + 



/^^ 



/ 



^u\{g^^-^a^'Pi{k-ujn)p,)u, 
2m 



x5{{k[ - fci + cjTi - ojt'Y - i?)e{uj\k[ - k^)) 



dri 



Ti — le 



(16) 



k is the meson momentum. The coupling constant g stands for the vector vertex wy'^u, whereas 
/ corresponds to the derivative coupling. The vector meson exchange generates its own vector 



contact interaction, corresponding to the cross on the vector meson hne. This is an additional 
contact interaction in the vertex NNjV. As indicated in 0, the contact term can be taken into 
account by the replacement, in the numerator, of the momentum k by the difference k — ujt, 
where ujt is the momentum of the spurion line connecting the ends of the vector meson line. 
This replacement has been done in (|T^). The factor ujti everywhere in {k — ujti) incorporates 
the contact terms both for the vector meson exchange and for the derivative coupling. We 
express k in terms of the lower vertex momenta: 

k — uoTi = ki — k[ — LOT for uj-{ki — k'^^) > 

k — uJTi = k[ — ki — ujt' for uj-{k'i — /ci) > (17) 



The first line of (|T^ corresponds to the first item in (|T6[), the second line corresponds to the 
second one. 

By this way, we find the vector contact term (both the left and the right one). It is 
convenient to represent it as the sum of two contributions: 

jur'"' = j'i^u' + j^y ■ (18) 

The first one, J^J^ ^ arises from the contraction of the A^A^-meson vertices with —gap 
in the meson propagator appearing in (0). For example, Jieft is obtained by the following 
substitution in eq. ([TlD (and similarly for Jright)'- 

V = g^-- J-a''\ik, - k[ - ujcpiT, r'))«', (19) 

V -> g^a + 7r-0^/3'a^(A;i - k[ - LUifiT, T')f. (20) 

2m 



The vertex V^ is given by ( [20| ) with a minus sign; this sign is compensated by the minus 
from —gap, what gives the substitution ([20|). In eqs. (p!9| , pOD we introduce the function, which 
takes into account the condition (|T^: 

^{t,t') = 1 ''; 'I ''^''\ (21) 

^^ ' ^ [ -r', \i X <x' ^ ' 

Since the vertex V^ appears multiplied on the left by {k\ —m,) and on the right by {k[ —m,), 
it can be transformed into: 

V'-^{g + fha - ^(^1 + K)a - i^ap,aUJ^'^{T, t') (22) 

The second item J^Jp in ([T8|) arises from the contraction of the vertices with the factor 
{k — ujTi)a{k — ujTijp in (|T^. It is given by (jT^) with the following vertices 

V^^{ki-k\), V ^ -uj^{t,t'). (23) 

/x /i 

The other terms, like the term proportional to cj, do not contribute to the upper vertex, 
since being multiplied by Co from the contact term it gives zero. The difference ki — k[ in the 
lower vertex appears in the expression {ki — m){ki — k[){k[ — m) = and, hence, does not 
contribute as well. The minus sign from Uc'yiUc = —7/3 is incorporated in the function (|2l|) . 



5 Calculating form factors 



The next step consists in extracting the deuteron form factors from the electromagnetic vertex 
Jjl^. As already mentioned, in contrast to the wave function which is always off energy shell, 
the on-shell amplitudes should not depend on the light-front plane orientation. However, in 
practice, due to the incompatibility of the transformation properties of the approximate current 
and wave function, the nonphysical u dependence survives in the on-energy shell deuteron 
electromagnetic vertex. Thanks to covariance, the general form of this dependence can be 
found explicitly |19[ : 
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^l9fj.iy + ^S 
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2uj-p 



2uj-p 



+ 



2uj-p 



(24) 



Here q = p' — p, P = p + p', with p and p' the initial and final deuteron momenta. As an 
example, we keep in (^4]) three tu-dependent terms only, though the total number of them is 
eight. Since we assume uj-q = (what corresponds to g+ = in the standard approach), all the 
form factors in (^4]) depend on Q^ = —q^ only. 

The a;-dependence of the wave function is not the only source of the nonphysical con- 
tributions -Bi-s- With an n- independent wave function (neglecting /s in (Q)) we still get 
the a;-dependent deuteron electromagnetic vertex. On the contrary, the physical form fac- 
tors J^iiJ^2iQi do not depend on uj even with an cj-dependent wave function. We emphasize 
that the physical tu-dependent extra components in the wave function (^) and the nonphysi- 
cal cj-depending structures in (^) are present not because of the covariant formulation of the 
LFD. Their counterparts appear in the non-covariant light-front approach too but the covariant 
approach allows to indicate them explicitly. 

The explicit formulas to extract the physical form factors Ti, T21 Q\ from the electromag- 
netic vertex and to separate them from nonphysical contributions -Bi_8 were derived in refs. 

As an example, we give here the formula for T\ only: 



T^ 



jp ^P 

'"'2UJ-P 



fj,iy 



q^'q" 



2uj-p 



p,,," 



P' 



uo'^uo 



A{uj-py' 



(25) 



The contraction with ujp in (^) corresponds in the standard approach to the component J^^. 
The expression for JF2 is proportional to the same contraction. However the expression for Qi 
contains not only the contraction with uOp but also others terms. Hence, the J^ component is 
not enough to find Qi. 

One can easily check, using uj-q = 0, that one has for the nucleon electromagnetic vertex 
(I) il^ujpTf = UpTPcj = 0. Here u appears from the contact interaction, the contraction UpT^ 
is the result of eq.(p5|). Therefore the contact terms do not contribute to Tx.T^-, due to the 
contraction with LOp in (pS]). They contribute only to Q^. The form factors found this way are 
still not exact, but they are separated from the non-physical ones -Bi-s- As shown in |jl9| and 
in [^, this separation is important, since the admixture of Bi can considerably change the 
results. 

After calculating traces and contracting tensors in eq. ( PSJ ) for T\ (and similar expressions 
for T2 and for Q^) we obtain the integrand depending on the scalar products of the four- 
vectors C(J, fci, k^i k[, k2, p, p', q with each other. In their turn, these scalar products, as well as 



the arguments of the wave functions, are expressed through the integration momenta and the 
momentum transfer. Their exphcit kinematical formulas are given in [^]. The integrands also 
contain the scalar products with uj but after integration this dependence disappears. 

By this way, substituting J^^^ in eq.(p5D (and in similar expressions |jl9|] for ^^2,^1)5 we 



find the analytical expressions of the form factors integrands. Then we carry out the numerical 
integration which, for the contact terms (two-loop), is 6-dimensional. 

6 Influence of the nucleon EM form factors 

Together with the deuteron wave function, an essential ingredient in these calculations are the 
nucleon electromagnetic form factors (NEMFF). The results presented in [Q were obtained 



with the dipole parametrization taken from Bilenkaya et al. |2T[ denoted hereafter by BKL. 
We analyse in what follows the influence of different NEMFF parametrizations in calculating 
the deuteron structure function within the same theoretical scheme. The parametrizations we 



consider, used in similar calculations found in the literature [0,[T3, 12^-03, are those given by 



Galster et al. ||2^, Hohler et al. p3[, a combination of Simon et al. for proton and Platchkov 
et al. for neutron (SP) |2^ and the more recent publication of Mergell et al. (MMD) [p5|] . 
The comparison [^ of the proton (G^;, G\j) and neutron {0%, G\[) charge and magnetic form 
factors with the existing experimental data [^ is shown in figure |^ . The values of the neutron 



charge form factor have been squared to account for the results reported in [30|. A first sight to 
this figure shows that some of these form factors have been apparently determined in the low 
momentum region and can not be used above 1-2 (GeV/c)^. The more adequated parametriza- 



tion covering the whole momentum region seems to be MMD ||25|, which will be adopted in 
our calculations. It provides an acceptable description for all the form factors although some 
variations can not be excluded due to the inaccuracy of the existing measurements, especially 
for the neutron form factor (see below). Any comparison between the theoretical and exper- 
imental results has to be understood in the context of this uncertainty that we would like to 
estimate. 

The influence of different NEMFF parametrizations on the deuteron structure function 
A{Q'^) is shown in figure |[ The results in figure H(a) are fully non relativistic whereas those 
in ^(b) correspond to LED impulse approximation. The choice of the NEMFF has small 
infiuence at low momentum transfer. For instance at Q^ = 0.5 (GeV/c)^, the deviation in the 
nonrelativistic structure function A is ~ 15%. This deviation increases with Q^ and reaches 
already a factor two at Q"^ = 2 (GeV/c)^ between BKL and SP parametrizations. The difference 
becomes dramatic for the last measured values Q^ = 6 (GeV/c)^ where there is one order of 
magnitude in the structure functions ^(Q^) calculated using BKL or Hohler parametrizations. 
It is worth noticing however that the Galster and MMD sets lead to close results in all the 
momentum range. Their difference in the LFD-IA relativistic calculation of A remains at the 
level of 10 % for Q^ < 10 (GeV/c)^. Similar deviations have been found in the structure 
function B. 

As one can see from figure |^, the neutron electric form factor G^ contains the greatest 
uncertainty. To estimate the infiuence of this uncertainty, we have considered for Q^ = 2 — 4 
(GeV/c)^ the maximal G^^ values compatible with the error bars, which are a priori not 
excluded by the experimental data although they considerably differ from all fits. The values 
of AIQ"^) obtained with such a modification exceed by 20-30% those obtained with the original 
MMD parametrization. From these considerations it follows that i) the comparison between 



Nucleons form factors, param SPA 



S^ 1.4 
1.2 
1 
0.8 
0.6 
0.4 
0.2 




- 






,,•"■* 




- 




,•*" 


■■■■ 






- 


.-•••:^ 


'^ 


r ^ 


pF"«Sr ■ 





* ■ > • . 


• -4?. .T 


rJ 


^ 


















: / 










:/ 










/ 


1 1 1 


1 1 1 


1 1 1 


1 1 1 



o 


0.5 
0.45 




0.4 




0.35 




0.3 




0.25 




0.2 




0.15 




0.1 




0.05 








E 


..-•• 




■' / 


.• 




E 




^'"'P'^ 






-- 




E 


••* 


v 


1 


» 


I 


■•i-.-f t ^: 


y 








E /,, 


k 


^-' 










•® 


r^ . 


r 










j^ 


























eT 












f 1 1 1 1 1 1 1 


1 1 1 1 1 





Q^ (GeV/c) 



10 

2 



Q (GeV/c)" 



o 



Coward (1968) 
Bartel(1973) 
Hanson (1973) 
Hohler (1976) 
Rock (1982) 



Q^ (GeV/c) 



10 

2 




Bosted(1990) 
A Markowitz (1993) 

* Lung (1993) 
■ Sill (1993) 

* Nikhef93 



6 8 10 

Q^ (GeV/c)^ 

+ Meyerhoff (1994) 

• Handivahis (1994) 

* Bruins (1995) 



Figure 3: Comparison between proton {G\ j^) and neutron {G% jyj) electromagnetic form fac- 
tors parametrizations and experimental data taken from [^ . The continuous line corresponds 

to MMD [|5l, dot-dashed (— -) to Hohler et al. Hi, dashed ( ) to Galster et al. [||; 

short-dotted (...) to BKL El and long-dotted (. . .) to SP 
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Figure 4: Influence of the nucleon electromagnetic form factors in the deuteron structure func- 
tion A(Q^) in the non relativistic impulse approximation (a) and in the LFD impulse approxi- 
mation (b) 

different calculations is only meaningful when the same parametrization of NEMFF is used and 
ii) the comparison between a theoretical prediction and the measured deuteron form factors 
is limited by the considerable uncertainties implied by the poor knowledge of the NEMFF for 

g2 > 1 (GeV/c)2. 



7 Results and discussion 

We present in this section the comparison of the deuteron structure functions and tensor polar- 
ization observable with the last measured values. In terms of the deuteron form factors these 
observables read 

2 
-? 

4 

^\i]FcFq + %fFl + ^^ (l + 2(1 + ri) tan^ U ^ -^ 



t20 [A{q^) 
where 



A{i 



B( 



F'c{^') + h'FQ(^') + hFl'"' 



9 



tan=^ -e B{q^\ 



^ 



(26) 
ri{l+v)FW). (27) 

t^FcFq + Iri^Fl + ^r^ ( 1 + 2(1 + 77) tan^ ^-9 ] F'i, 



{2i 



'c 



2rj 



-^1- ^[^1 + ^1 -^2(1 + ^7)] 



Fm 
Fn 



-^1-6^1+^2(1 + ^) 



and 7] = Q'^/AM'^. In comparing our calculations with the data we should first emphasize that, 
though the theoretical framework we use - the light-front dynamics - is fully selfconsistent, 
the validity of our perturbative method to evaluate the wave function is restricted to relative 
nucleon momenta k smaller than the nucleon mass m. Since in a form factor calculation the 
momentum transfer is distributed between two deuteron vertices, the initial and final ones, this 
value of k corresponds approximately to the momentum transfer Q"^ < (2r7i)^ f» 3.5 (GeV/c)^. 
In absence of subtle cancellations increasing the uncertainty in the theoretical predictions, we 
can expect a reasonable description of the experimental data in this momentum region. 
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Figure 5: The structure function A(Q^) of the deuteron with MMD nucleon form factors. The 
curves are explained in the text. 

In figure ^ the structure function v4((5^) is shown together with the experimental points 
recently obtained at TJNAF (Hall C in filled circles |^ and Hall A in opaque triangles up |0]). 
The dashed curve corresponds to the non-relativistic impulse approximation with the S- and D- 
waves of the Bonn-QA wave function 0. The dot-dashed line is also calculated in the impulse 
approximation but using the light-front formalism with relativistic deuteron components /i and 
/2 and the /s component in first degree only. The solid line incorporates, in addition, the contact 
terms which turn to have small influence in this observable. We can see a good agreement until 
Q^ = 2.5 (GeV/c)^, a momentum region where the departure from a non relativistic description 
reaches one order of magnitude. It is worth noticing that this agreement is found with the LED 
wave function alone, i.e. without explicitly including any MEG diagrams. However, as it has 
been already mentioned, the extra component /s accounts for the so called pair terms in the 
deuteron electrodisintegration amplitude JlSl . We remark that a systematic deviation seems to 
manifest above 3 (GeV/c)^. This excess is smoothly increasing with increase of Q^. However, 
this is the region where the perturbative calculation is hardly reliable, since the deuteron wave 



function appears in v4(Q^) in the fourth degree and a 50% correction may change the value by 
a factor 5. Besides, as it has been shown above, at Q^ > 2 (GeV/c)^ the uncertainty coming 
from the NEMFF is also high. 
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Figure 6: The deuteron tensor polarization ^20 at 9 
the same as in figs. ^. 



70°. The designations of the curves are 



The same agreement is seen in the deuteron polarization observable ^20 displayed in figure 
^ using the same drawing conventions. Again a sizeable effect of the relativistic corrections 
is properly taken into account in our calculations although the actual error bars are quite 
comfortable. 

The calculated structure function B{Q^) is shown in figure |^, using the same drawing 
conventions, together with the experimental data from |3ll]. Comparing the B[Q^) results with 
those of v4(Q^) and ^20? one can see a considerable deviation of our B calculation from the 
data, especially in the region of the minimum, as well as highest sensitivity to the different 
approximations. This minimum corresponds to the zero of the deuteron magnetic form factor 
Fm, and a small shift on this value drastically changes the value of B{Q'^) in its neighborhood. 
It is important to emphasize that the zero of Fm which exists in a nonrelativistic calculation, 
disappears for the relativistic /i, /2 and appears again when /s is taken into account p. In a 
relativistic framework this minimum is thus a consequence of a delicate cancellation between 
the /s contribution with the contributions of /i and /2. None of these /i,/2,/5 components 
has been calculated with enough accuracy in our perturbative approach. Therefore we cannot 
pretend to a detailed description of B^Q"^) in this region. 

The situation is different for AiQ"^) and ^20- Figure [^ shows the relative contributions of 
the three deuteron form factors to these observables according to equations ( p6|) and (|^). One 
can see that their values are dominated by the charge and quadrupole form factors Fc and Fq, 
whereas the contribution of Fm, containing the highest uncertainty, is suppressed. 

Taking into account the same order of accuracy, i.e., f?^, which was kept in the calculation 
of the wave function, we get a good description of the data on ^20 as well as of A[Q'^) in the 
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Figure 7: Same as figure ^ but for -B(Q^ 



same region of tlie momentum transfer, up to 2 (GeV/c)^. Tfie relativistic effects in /i, /2 and 
tlie contribution of /s are important to achieve tliis good description. Tlie contribution of tfie 
contact term to A{Q^) and ^20 is small, since, as noted above, it contributes to Fm only, whereas 
for B it is small in the momentum transfer region shown in figure |^ but increases at higher 
momenta (e.g. by a factor 2 at Q^ = 6 (GeV/c)^ ). These results show that the relativistic 
effects in the deuteron wave function, including the extra component /s, and in the deuteron 
electromagnetic vertex make considerable influence on the deuteron form factors. Note that 
the effect of incorporating /s on ^20 is qualitatively similar to that obtained when adding the 
contribution of the pair current in non-relativistic calculations P? . 



We should finally emphasize that the results presented in this section have been obtained 
without fitting any new set of parameters, neither in the nucleon-meson form factors nor in the 
NN interaction kernel, which was taken as it is given in ref. 0. 

The agreement of our calculations with the experimental data, as well as with some 
other calculations carried out in the framework of meson- nucleon dynamics (see e.g. [[7|-p7||), 
shows that the deuteron structure at small distances is understood rather well within this 
theoretical framework. It is a remarkable fact that probing distances of the order of 0.1 fm, at 
which quarks effects should manifest themselves in their full glory, could be accounted by the 
relativistic nucleon-meson dynamics and the phenomenological nucleon form factors. 

The parameters of the effective nucleon-meson Lagrangian, fixed by fitting the A^A^ ex- 
perimental data, should be in principle derived from QCD, what is a separate problem. This 
allows us to make the following general conclusion: at least in what concerns the deuteron, the 
relativistic nuclear dynamics can be developed independently of its derivation from QCD. 

Our calculations can be improved in many respects. An exact solution of the equation for 
the deuteron is necessary and is in progress. Besides the non-perturbative calculation, it implies 
the determination of the A^A^ kernel parameters in the framework of the light-front equations. 
This requires also a careful treatment of higher order contributions to the kernel. Finally, the 
calculation of electromagnetic observables should include higher Fock states {NNtt, etc.) and 
the meson exchange currents which are not included in the wave function components, like 
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Figure 8: The contributions of the different deuteron form factors {Fc-, Fq, Fm) to the structure 
function A(Q^) (a) and ^20 (b). 

those corresponding to the interaction of the photon with the intermediate mesons (p7r7, . . .). 
The improvement of the experimental data on the nucleon EM form factors is also an 
urgent task. 
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